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ABSTRACT 

We consider the excitation of the inertial modes of a uniformly rotating fully convective 
body due to a close encounter with another object. This could lead to a tidal capture 
or orbital circularisation depending on whether the initial orbit is unbound or highly 
eccentric. We develop a general self-adjoint formalism for the response problem and 
thus solve it taking into account the inertial modes with m = 2 for a full polytrope 
with n = 1.5. We are accordingly able to show in this case that the excitation of 
inertial modes dominates the response for large impact parameters and thus cannot 
be neglected in calculations of tidal energy and angular momentum exchange or orbital 
circularisation from large eccentricity. 
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1 INTRODUCTION 

The process of tidal capture of a body into a bound or- 
bit through the excitation of oscillatory internal modes is 
thought to be of general importance in astrophysics. It is 
likely to play a role in binary formation in globular clusters 
(eg. Press & Teukolsky, 1977, hereafter PT and references 
therein) as well as in tidal interactions of stars within galac- 
tic centres. When there are repeated close encounters as for 
a highly eccentric orbit, circularisation may result, a process 
of potential importance for extrasolar planets (eg. Ivanov & 
Papaloizou 2004, hereafter IP). 

Until now, only the oscillation modes associated with 
spherical non rotating stars have been considered with pos- 
sibly rotation being treated as a perturbation (eg. PT,IP). 
But angular momentum is transfered in tidal encounters and 
when, as is often the case, the internal inertia is much less 
than that of the orbit, the tidally excited object would be 
expected to rotate at a balanced rate or undergo pseudo syn- 
chronisation. Then the rotation frequency is matched to the 
most important ones in a Fourier decomposition of the tidal 
forcing. It is then natural to expect that inertial modes con- 
trolled by the rotation (see eg. Papaloizou & Pringle 1981, 
hereafter PP) become important for the tidal response. This 
is especially the case for a barotropic configuration 1 , with 
no stratification and hence no low frequency g modes. We 



* E-mail: J.C.B.Papaloizou@damtp.cam.ac.uk (JCBP); 
P.Ivanov@damtp.cam.ac.uk (PBI) 

For a barotropic fluid the pressure P has the same functional 



would expect inertial modes, which have periods that scale 
with the rotation frequency, to be particularly important 
for large periastron distances as the fundamental spherical 
modes suffer increasing frequency mismatch as this distance 
increases. 

It is the purpose of this paper to calculate the response 
of a rotating body taking the inertial modes into account. In 
particular we focus on modes with m — 2 which dominate 
the angular momentum exchange and a full polytrope with 
index n = 1.5. We show that the inertial modes indeed dom- 
inate the response at larger periastron distances of interest 
in many cases. 

In section 2 we give the basic equations governing the 
linear response of a uniformly rotating barotropic star to 
the tidal forcing due to a close encounter. In section 3 we 
give a self-adjoint representation of the linear problem and 
discuss how it can be used to provide a formal solution in 
terms of a spectral decomposition as can be done for a non 
rotating star. In section 4 we develop expressions for the 
energy and angular momentum exchange. In section 5 we 
give results for a polytrope with index n = 1.5, showing 
that the m — 2 inertial wave response is contained within 
just a few global modes and that inertial modes dominate 
the tidal response for large periastron passages. Accordingly 
these must be considered when discussing tidal circularisa- 
tion from large eccentricity for objects such as extrasolar 
planets. 

dependence on the density p: P = P(p), for the unperturbed 
configuration and during perturbation. 
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2 BASIC EQUATIONS 

We assume that the star is rotating with uniform angular ve- 
locity ft directed perpendicular to the orbital plane. In this 
case the hydrodynamic equations for the perturbed quan- 
tities take the simplest form in the rotating frame. For a 
convective star, approximated as barotropic, we have 



Dt 2 Dt 



-VW. (1) 
Here £ is the Lagrangian displacement vector, and 



W = c s p /p + * lnt +* e 



(2) 



where p is the density, p is the density perturbation, c s is 
the adiabatic sound speed, v]/ lnt is the stellar gravitational 
potential arising from the perturbations and ty ext is the ex- 
ternal forcing tidal potential. Note that the centrifugal term 
is absent in equation Q being formally incorporated into 
the potential governing the static equilibrium of the unper- 
turbed star. The convective derivative -St = Mr as there is 

Dt at 

no unperturbed motion in the rotating frame. 

The linearised continuity and Poisson equations are 



p = — V(p£) and 



A*" 1 * = 4nGp 



(3) 



We use a cylindrical coordinate system (uj, <f>, z) with 
origin at the centre of the star and assume that 
the tidal potential and the resulting perturbation to 
some quantity, say Q, is represented in terms of a 
Fourier transform over the angle (p and time t. Thus 

Q — y"] m (^exp(im(f)) f + .°° dcyQ m exp(— iat) + ccj , where 

the sum is over m = and 2 and cc denotes the complex con- 
jugate. The reality of Q implies that the Fourier transform, 
indicated by tilde satisfies Q m (<j) = Ql m (— &). 
To solve the forcing problem we obtain equations for the 
Fourier transforms of the perturbations using Q — We 
can then express £ m in terms of W m with help of Q, the 
density perturbation in terms of and W m with 

help of and substitute the results into the continuity 
equation to get 

a 2 AW m - oBW m - CW m = a 2 d4(*- ' + - W m ), (4) 



where d = 4Q. 2 — a 2 , and 











dz 





+ 



m 2 p 



■uj dvj dz 



(4) 



(5) 
(6) 



It is very important to note that the operators A, B and C 
are self-adjoint when the inner product 



(Wi\W 2 ) 



dz-ujd-mWiW2, 



(7) 



with V denoting the volume of the star. Also when W is not 
a constant A and B are positive definite and C is non neg- 
ative definite. Equations © anQl ^ ne Poisson equation 
$Q form a complete set. 

When = 0, equation (0J describes the free oscillations 
of a rotating star. These may be classified as relatively high 
frequency / and p modes and inertial modes with eigen fre- 
quencies a ~ fl. The / and p modes exist in non rotating 



stars and can be treated in a framework of perturbation the- 
ory (eg. IP) . Here we focus on the response of inertial waves 
to tidal forcing, for slow rotation. We seek a solution as a 
series expansion in the small parameters a/Q t and f2/f2», 
where Q t = */ (GM*/R%), with M* and R t being the mass 
and the radius of the star, respectively. The small parame- 
ters will be assumed to be of the same order. 
To zeroth order (when a — > 0) the tidal potential induces a 
quasi-static bulge in which tidal forces are balanced by pres- 
sure and self gravity and p m — > p m ,st with 5'™' — > &m%t- 
From equation (|1J we conclude that W is smaller than the 
potential perturbations by a factor of order a 2 /Q. 2 (see PP 
for details). Accordingly, setting W — in equation @ and 
using equation (J^J we obtain 

P 



Pm.st 



■ \ * r. 



A* 



4:TTGp m ,S 



(8) 



In the lowest order approximation we find Wm using equa- 
tion |Q} neglecting the higher order term proportional to 
pW m /c 2 on the right hand side. We thus have 



a 2 AW m - aBWm - CW m = S, 

where the source term 

S = — o 2 dp m , st = a 2 d^(^'^ t + SS?l 



(9) 



(10) 



is completely determined by the external forcing via equa- 
tion (JHJ. Finally, the perturbation of the internal gravita- 
tional potential associated with W m , 59 m = — ^>m% 
follows from © and JHJ 



A<5* r 



(11) 



3 REDUCTION OF THE PROBLEM TO 

STANDARD SELF-ADJOINT FORM AND 
FORMAL SOLUTION 

Equation © is of a general type that governs the linear re- 
sponse of rotating bodies in a context that is much wider 
than the specific one considered here ( see Lynden Bell & 
Ostriker 1967). One recurring difficulty has been that the 
way in which the frequency a occurs does not readily allow 
application of the general spectral theory of self-adjoint op- 
erators to the response calculation even though A, B and 
C are themselves self-adjoint (eg. Dyson & Schutz 1979). 
We here provide a formulation in which the problem can be 
dealt with in the same way as a standard self-adjoint one 
such as occurs for a non rotating star. 

We begin by remarking that the self-adjoint and non nega- 
tive character of A, B and C allows the spectral theorem 
to be used to specify their square roots, eg. A 1 / 2 , defined 
by condition A = A 1 ^ 2 A 1 * 12 . The positiveness of A allows 
definition of the inverse of A 1 ^ 2 , A~ 1//2 . Technically, this is 
done through using the spectral decomposition 



AW = 



i 



(W\*i) F (A)W = F (^Hj {Wfa) , (12) 



..etc., where 4>i an d Aj are the eigenfunctions and the eigen- 
values of the corresponding operators and the last equation 
defines a general function of the operator A. The sum may 
formally include integration with an appropriate measure 
in the case of a continuous spectrum of the corresponding 
operators. 
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Now let us consider a generalised two dimensional vector 
Z with components such that Z = (Zi,C 1 ' 2 W m ) and the 
straightforward generalisation of the inner product given by 
equation It is easy to see that equation © can be de- 
rived from the matrix form 



aZ = HZ + S, 
where 

A-^BA- 1 ' 2 



H 



c l/2 A -l/2 



A -l/2 C l/2 




(13) 



(14) 



and the source vector S has the components (A^ 1 ^ 2 S, 0). It 
then follows from (1131 that we may choose Z\ = a A 1 ^ 2 W m - 
Since the off diagonal elements in the matrix (1141 are adjoint 
of each other and the diagonal elements are self adjoint, that 
the operator H is self-adjoint is manifest. Now we can use 
the spectral theory of self-adjoint operators to look for a 
solution to J3J in the form Z — ctkZk, where Zk are the 
real eigenfunctions of H which satisfy 



OkZk — HZk 



(15) 



the associated necessarily real eigen frequencies being Ok- 
Substituting 1151 into 1131 we obtain 

< z k \s > 



ctk = — - — , (16) 

< Zk\Zk > {o + iv — <Tfc) 

where, following the Landau prescription, we have added an 
infinitesimal imaginary part v > to the frequency a. The 
inner product induced by H turns out to be < Zk\Zi >= 
a k ai(W h \AW,) + (W k \CW,), where W k = C~ 1/2 Z^ = 
a~ 1 A~ 1/2 Z'[. Using (EHJ and JTSJl we explicitly obtain 

W m = adY j f) {W ^ M . W k , (17) 
*-~t N k (a k -IV - a) 

k ' 

where N k = <j 2 k {W k \AW k ) + {W k \CW k ) is the norm. 

4 TRANSFER OF ENERGY AND ANGULAR 
MOMENTUM FROM ORBIT TO INERTIAL 
WAVES 

The time derivative E of the wave energy E associated with 
the rotating frame and the time derivative L of the wave 
angular momentum L are given by usual expressions 



E 



d(j> ( p |>]> e:!;t 



' d<j> 



(18) 



The wave energy associated with the inertial frame, Ei is 
given in terms of E and L (eg. Friedman & Schutz 1978, 
hereafter FS) through 



Ei = E + QL. 



(19) 



Expressing p with help of © in terms of W and dis- 
carding the total derivatives which cannot contribute to the 
total energy transfer AE = J_ dtE, we obtain 



AE 



dt / , 

•oo JO 



(4(w"-^)i* e 



(20) 



Contributions AE m to AE arising from different values of m 
may be calculated separately. Substituting the correspond- 
ing Fourier transforms into 1201 . and integrating over time, 
we obtain 



AE„ 



47r 2 i 



daa{{^T\5p m ) ~ («5p m |*r)}, 



(21) 



where 8p m — -%z(W m — S4> m ). Note that for the problem on 
hand an d are rea l ( see below). Now we substi- 

tute the spectral decomposition <5^/ m = ^2i k &k5^k, where 
A5* fe = ^£-{Wk-8<l!k), together with JTTJ into (STJ. The 
singularities in the integral over a corresponding to a — Ok 
and determining all non-zero contributions to the integral 
are dealt with using the Landau prescription. Then the en- 
ergy transfer is a sum of contributions associated with each 
mode, being given by 

AE m = -8^ 3 V g fc( 4 " 2 ~ (W k \p m ,st)(6 Pm . k \*T), (22) 
* — ' N k 

k 

where Sp m ,k = -4" (Wit — ^k)i and both integrals are evalu- 
ated for a — Ok- 

Using the symmetry properties of the operator A + 4 ^ p one 
can readily see that the integrals are, in fact, the negatives 
of each other, thus we finally obtain 



AEr, 



8tt 



3V a U^ 2 -<7k) ,- 



Si 



N k 



At, 



(23) 



where A k = \{p m , s t\Wh)\- 

The expression for the angular momentum transfer follows 
immediately from (1231 and the general relation between 
wave energy and angular momentum corresponding to a 
mode with particular value of m: E m /L m = er/m (eg. FS). 
We have ALo = and 



AL 2 



16tt 



a 3 k (in 2 



k 



N k 



-Al 



(24) 



The expressions 1231 . 1241 can be evaluated when the 
star is in a parabolic orbit, being the limit of a highly ec- 
centric orbit, using the expression for ^i ext given in the 
quadrupole approximation by PT. From this we obtain 



5 (2-|m|)! B ]ml 
2tt 3 (2 + jmj)! (1 + 9) 



h,- m (y)n p r 2 Pl ml (cos 9), (25) 



where B 2 = Vlf > B° = — Ivf " q is tne ratio of M * to 
the mass of the star exerting the tides, M: q = M*/M. 



GM(l+i) 



is a typical frequency of periastron passage 

and R p is the periastron distance. The functions l2,-m{y) 
specify the dependence of 9%%* on a, they are described by 
PT. y—fl(a + m), where we use dimensionless frequencies 
a — a/Q and Q — Q/Q p . P™(x) is the associated Legendre 
function, and (r, 9) are the spherical coordinates. 
Now we express the density perturbation entering in 1231 . 
1241 as p~m,st = — -^j^^m • Here we introduce a factor F 
to account for self-gravity of the star that can be obtained 
from the solution of @ 2 . We substitute this expression 
and J2U in gSJ, CU and express all quantities in natural 
units. These are such that the spatial coordinates, density 
and sound speed are expressed in units of R t , the averaged 

3M, 



density p ■■ 



and 



Rt 



, respectively. We then have 



F = 1 in the frequently used Cowling approximation. 
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Figure 1. The dependence of on the values of o^. 



AL 2 



2C 2 



{^(4-^)Ql/ 2 2 _ m (y)}^, (26) 



« 3 X^( 4 -^)^4-m(y)}^r> (27) 



(i + g) 2 

where E» = L* = M*y/GM*R*, and 77 = fi*/fi p = 

(l+q)UB? ' ^ = B and = 3 / 4 ' Tne over l a P integrals 
Qfc have the form: 

Q fc = (4 i?r2 - p 2 m |^)/v / ^, (28) 
where jV fe = iVfe/fi 2 . 

5 NUMERICAL CALCULATIONS 

The eigen frequencies <Tk and the integrals Qk in equations 
126I27H have to be calculated numerically, for a given model 
of the star. Therefore, equation l|15p describing free oscilla- 
tions of the star must be solved. For our numerical work we 
model the star as n = 1.5 polytrope and use the method of 
solution of l|15|l proposed by PP 3 . Namely, we reduce l|15> to 
an ordinary matrix equation using an appropriate set of trial 
functions, and solve the corresponding eigenvalue problem. 
The details of our approach will be given in a separate publi- 
cation. Here we note that we use different numbers N t of the 
trial functions to study convergence of our numerical scheme 
with the largest number being equal to N™ ax = 225. Also we 
adopt a spherical model of the star assuming that the distor- 
tion of the stellar structure due to rotation does not influence 
significantly our results. Since the modes corresponding to 
m = 2 fully determine the transfer of angular momentum 
in the quadrupole approximation and apparently dominate 
the energy transfer, which can only be increased if m = 
modes are included, we consider only them. 

The absolute values of the eigen frequencies at are al- 
ways smaller than or equal to 2Q. For complete polytropes it 

3 In order to check the method we also calculate the eigen spec- 
trum of an n = 1 polytrope. Our results are in a good agreement 
with results obtained by Lockitch & Friedman (1999) and Din- 
trans & Ouyed (2001) who used different methods. 





A, 



Figure 2. The distribution of Wk over the star. W& = at the 
origin. Contour lines of different style correspond to values of 
Wk of opposite signs. The upper left (right) plot corresponds to 
ai = —1.063 ( <T2 = 0.697). The lower left (right) plot corresponds 
to 0-3 = 0.362 (<74 = 0.703). 




Figure 3. The dependence of the dimensionless quantities e, ei 
and A on Q (see 1291 and the text for definitions). 

is expected that they ultimately yield an everywhere dense 
discrete spectrum (PP). In our numerical approach their 
number is equal to 2Nt, and therefore the number of the 
terms in the series in 126H and 127H depends on Nt- How- 
ever, only the eigenfunctions with some definite values of 
o"fe give a significant contribution to the series. These terms 
correspond to global eigen modes with a large scale distri- 
bution of Wk over the star. Therefore, the overlap integrals 
corresponding to the global modes are not averaged to small 
values after integration over the volume of the star. The 
number of the global modes does not depend on Nt. In fig- 
ure 1 we show the overlap integrals as functions of the eigen 
frequencies. There is a sharp rise of Qk near <7fc = and 
very close to at = 2, and several peaks corresponding to 
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Figure 4. The dependence of the energy transfer on rj. The black 
curve (red and green) give the contribution of the inertial waves 
(fundamental mode). We assume that ij<1 and f2 = (l cr it ~ 1.6. 

different values of The rise at ou ~ 0, 2 is due to numer- 
ical inaccuracy of our method, but the corresponding modes 
practically do not influence our results and therefore the nu- 
merical inaccuracy is not significant for our purposes. Note 
that the feature at <Tk « can be removed by an appropriate 
regularisation of C. 

Three peaks with o\ = —1.063 and Q\ — 0.037, 
<T2 = 0.697 and Q 2 = 0.175, cr 3 = 0.362 and Q 3 = 0.034 
correspond to the global modes. The spatial distribution of 
the corresponding eigenfunctions is shown in figure 2. For 
comparison we also show a non global mode with 04 = 0.703 
and a very small value of Q4. It is clear from this figure that 
the global modes have a large scale distribution of Wk- 

It follows from 12611 and I27H that for a given model of 
the star, the quantities 



(l + < 



frf{AE m /E, 



), A = (l + g )V(ALa/L.),(29) 



depend only on Q. In figure 3 we show the dependence of ti 
(the solid black curve) and A (the solid red curve) on Cl. The 
solid green curve shows the dependence of the energy trans- 
fer in the inertial frame, = £2 + OA. Note that the energy 
transfer in the inertial frame is negative when H > 2. The 
dotted curves show the dependencies of the corresponding 
quantities calculated with only two global modes, the ret- 
rograde mode with o\ = —1.063 and the prograde mode 
with 02 = 0.697. Figure 3 shows that the dotted and the 
solid curves are very close to each other. Therefore, only 
two modes are essential for our problem. The angular mo- 
mentum transfer is equal to zero when Q, — Q, cr it = 1-6, 
where e cr i t ~ 0.0065. The condition AL2 = may be eas- 
ily realised in an astrophysical situation where the inertial 
modes dominate the tidal response and the moment of in- 
ertia of the star is sufficiently small. In this case the system 
quickly relaxes to the so-called state of pseudo synchronisa- 
tion where AL 2 = 0. In figure 4 we show the dependence 
of AE/E, on rj for this case. The black curve shows the 
contribution of the inertial waves, 



AE 2 



0.0065 
(l + g)V 



E*. 



(30) 



The red curve shows the contribution of the / mode to 
the energy transfer calculated for a non-rotating star (PT) 
and the green curve shows the same quantity but calculated 
when Cl = ilcrit (IP) 4 - It is clear from this figure that when 
Tj > 8 the inertial waves dominate the tidal response. For a 
planet with Jupiter mass and radius circularising from large 
eccentricity to attain a final period of 6 days, 77 ~ 18 (IP), 
such that inertial modes dominate by a wide margin. 

6 CONCLUSIONS 

In this Paper we formulate a new self-adjoint approach to 
the problem of oscillations of uniformly rotating fully con- 
vective star and apply this approach to the tidal excitation 
of the inertial waves and tidal capture. The approach can be 
used in any other problem where oscillations with typical fre- 
quencies of the order of the angular velocity are important. 
It can also be extended to the case of convectively stable 
stars. The oscillations with higher frequencies can also be 
treated in framework of an extension of our formalism. 

We show that the tidal response due to inertial waves 
can be represented as a spectral decomposition over nor- 
mal modes. This formalism is general enough to be applied 
even when the spectrum is not regular and discrete. For full 
polytropes, only a few global modes with a large scale distri- 
bution of the perturbed quantities over the star contribute 
significantly to the response. We find general expressions for 
the energy and angular momentum transfer between the or- 
bit and inertial waves. It is shown that when the inertial 
waves dominate the tidal response and a state of pseudo 
synchronisation is reached the energy transfer is determined 
by simple equation 13011 . As follows from that equation the 
energy transfer does not decrease exponentially with the or- 
bital periastron and therefore at large periastra the energy 
transfer due to inertial waves dominates over the energy 
transfer due to excitation of the fundamental mode. This 
can be very important for the problem of tidal circularisa- 
tion of the extra solar planets (see above, IP and references 
therein), and can be of interest for other problems, such as 
eg. the problem of tidal capture of convective stars in glob- 
ular clusters. The application of our approach to extra solar 
planets will be discussed in a separate publication. 
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the conclusions of the paper. 



